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We have studied the B^ J/ipn^ and — > J/tpK^ decays using a 20.7 flj^^ data set collected 
with the BABAR detector. We observe a signal of 51 ± 10 B^ — > J/ipn^ events and determine 
the ratio ^ J/tpir^yBiB^ JpK^) to be [3.91 ± 0.78(stat.) ± 0.19(sj/st.)]%. The CP- 

violating charge asymmetries for the B J/ip tt* and B^ J/'iji K^ decays are determined to be 
= 0.01 ± 0.22(stat.) ± 0.01{syst.) and Ak = 0.003 ± OmO{stat.) ± 0.004{syst.). 

PACS numbers; 13.25.Hw, 13.25.-k, 14.40.Nd 
The decay — > J/'^tt^ is both Cabibbo-suppressed and color-suppressed. If the leading-order tree diagram 
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is the dominant contribution, its branching fraction is 
expected to be about 5% of the Cabibbo-allowed mode 
J/ipK^. A comparable prediction can be obtained 
with a simple model based on the factorization hypoth- 
esis Previous studies of this decay were performed 
by the CLEO @ and CDF [| collaborations. Signifi- 
cant interference terms between the suppressed tree and 
penguin amplitudes could produce a direct CP-violating 
charge asymmetry in the — > J/ij^ir^ decays at the 
few percent level Q. On the contrary, a negligible di- 
rect CP- violation is expected in the B^ J/ipK^ de- 
cays, because for h ccs transitions the Standard Model 
predicts that the leading- and higher-order diagrams are 
characterized by the same weak phase. 

In this paper we present a measurement of the ratio of 
branching fractions B{B^ J/ip7T^)/B{B^ J/ipK^) 
along with a search for direct CP-violation in these 
channels. The data were recorded at the T(4>S') res- 
onance in 1999-2000 with the BABAR detector at the 
PEP-II asymmetric-energy e~^e~ collider at the Stan- 
ford Linear Accelerator Center. The integrated lumi- 
nosity is 20.7 fb~^, corresponding to 22.7 million BB 
pairs. We fully reconstruct B^ J/iph^ decays, where 
= , . Signal yields and charge asymmetries 
are determined from an unbinned maximum likelihood 
fit that exploits the kinematics of the decay to identify 
the TT^, and background components in the sample. 
This kinematic separation is sufficiently good so that no 
explicit particle identification is required on the charged 
hadron h^, thereby simplifying the analysis. At the same 
time, particle identification can be used to perform a 
crosscheck of the measurement. 

The BABAR detector is described in detail elsewhere . 
A five-layer silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH), in a 1.5-T solenoidal magnetic 
field, provide detection of charged particles and measure- 
ment of their momenta. The transverse momentum res- 
olution is apjpt = (0.13 ± 0.01)% • pt + (0.45 ± 0.03)%, 
where pt is measured in GeV/c. Electrons are detected in 
a Csl electromagnetic calorimeter (EMC), while muons 
are identified in the magnetic flux return system (IFR), 
which is instrumented with multiple layers of resistive 
plate chambers. A ring-imaging Cherenkov detector 
(DIRC) with a quartz bar radiator provides charged par- 
ticle identification. 

An electron candidate is selected according to the ratio 
of the energy detected in the EMC to track momentum, 
the cluster shape in the EMC, the energy loss in the DCH, 
and the DIRC Cherenkov angle, if available. A muon can- 
didate is selected according to the difference between the 
expected and measured thickness of absorber traversed, 
the match of the hits in the IFR with the extrapolated 
track, the average and spread in the number of hits per 
IFR layer, and the energy detected in the EMC. 

J/^jj — > p.^ fJ.~ candidates are constructed from two 
identified muons with polar angle in the range [0.3, 2.7] 



radians and with invariant mass 3.06 < M^+^~ < 
3.14GeV/c^. The absolute value of the cosine of the 
helicity angle of the J/-0 decay is required to be less 
than 0.9. J/tp e+e" candidates are constructed from 
two identified electrons with polar angle in the range 
[0.41,2.409] radians and with invariant mass 2.95 < 
A/g+e- < 3.14GeV/c^. The absolute value of the cosine 
of the helicity angle is required to be less than 0.8. 
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FIG. 1: Distribution of AEk vs. AE.^ for 



J/xjjK^ and 



J/?/>7r* events from Monte Carlo simulations. 



B^ candidates are formed from the combination of 
a reconstructed J/'0 , constrained to the world average 
mass [P, and a charged track . A vertex constraint 
is applied to the reconstructed tracks before computing 
two kinematic quantities of the B^ candidate used to dis- 
criminate signal from background. We define the beam 
energy-substituted mass toes as 



"^ES 



[(./2 + Pi.pB)Vi^f]-|PB| 



(1) 



where ^/s is the total energy of the e+e~ system in the 
T{AS) rest frame, and (i?i,pi) and (-Eb,pb) are the 
four-momenta of the e+e" system and the reconstructed 
B candidate, both in the laboratory frame. We define 
the kinematic variable Ai5^ {/S.Ek) as the difference 
between the reconstructed energy of the B^ candidate 
and the beam energy in the T{AS) rest frame assuming 
/i± = TT±{K^). We require |A£:^| < 120 MeV, \/:^Ek\ < 
120 MeV and toes > 5.2 GeV/c^. Figure shows the dis- 
tribution for Monte Carlo simulations of B^ J/ipir^ 
and B"^ J/^K^ events in the {AE^, AEk) plane. 
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FIG. 2: The AEk distribution and fit for the events in the 



data sample with mES > 5.27GeV/c 
resents the background contribution. 



The dashed curve rep- 



The selected sample contains 1074 J/%p{-^ 
fi+fi-)h'^ and 1081 B± J/'0(-> e+e")/i=^ candidates. 
A fit to the AEpc distribution with the sum of a Gaussian 
and a polynomial function, modeling the — > J/^K^ 
signal and the background contribution, is shown in 
Fig. |. 
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FIG. 3: Distribution of AEk vs. AE-n for the events in the 
rriES sideband of the data sample. 



The background contaminating the sample is char- 
acterized with events in the data that are sufficiently 
far from the typical signal regions (sidebands of the 
data sample). We define ities sideband events by the 
requirement that 5.2 < mEs < Mb — 4cr(mEs) = 



5.27 GeV/c^ where Mb is the world average mass tJ 
and cr(rnEs) is the toes resolution; their distribution 
in the (AE.^, AEk) plane is shown in Fig. ||. We de- 
fine AEfc and AiJ^ sideband events by the require- 
ment that 120 > \AEk\ > ia{AE) = 42MeV and 
120 > \AE^\ > Aa{AE) = 42MeV, where a{AE) is the 
width of the fitted Gaussian in Fig. |^. The distribution 
in Toes of the sideband events is modeled by an ARGUS 
function , with an additional Gaussian peak in the toes 
signal region for events from other B — > J/ipX decays. 
The number of background events in this peak has been 
estimated to be 10 ±4 with detailed Monte Carlo simula- 
tion of inclusive charmonium decays. Figure ^ shows the 
Toes distribution for the data sample, along with the fit. 
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FIG. 4: The rriES distribution and fit for the events in the data 
sample. The ARGUS (dashed curve) and peaking (dotted 
curve) components of the background are also displayed. 



Our fit to the data sample is based on maximizing the 
following extended likelihood function: 



M 
j=l i 



mEs')N^, (2) 



where j is the index of the event, i is the index of the 
hypothesis {i = tt, K, bkd) , Ni are the yields for the 
B^ J/ipn^, B^ — > J/ipK^, and background events 
in the sample, and M is the total number of events. 
The observables AE^, the momentum p of the final-state 
charged hadron computed in the laboratory frame, and 
Toes are used as arguments of the probability density 
functions (PDF) Pi. The PDFs are mainly determined 
from data with limited input from simulation. 

It is useful to define the new variables D — AEk — 
AEt^ ~ 7 {^^p^ + m\ — ^'p^ + mj^ , where 7 is the 
Lorentz boost from the laboratory frame to the T(45) 
rest frame, and S = AEk + AE^, = 2AE.„ -I- D. These 
variables have the property that {AEt^, D) in the pion hy- 
pothesis, [AEktD) in the kaon hypothesis, and {S,D) 
in the background hypothesis are uncorrelated at the 
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1% level. Therefore, with appropriate transformations 
of variables, each Pi(Ai?7r,P, ^es) can be written as a 
product of one-dimensional PDFs: 

P^{AE^,p,mEs) = /7r(Ai;^)g^(D)/i^(mEs) , (3) 
PK{AE^,p,mEs) = fK{AEK)gK{D)hK{mEs),{4=) 
Pbkd{AET,,p,mEs) = fbkd{S)gbkd{D)hbkd{mEs)- (5) 

The f^{AE^), fK{AEK), K{m-Es), and /lif(mEs) 
components are the AE and tties resolution functions for 
the signals. The mean values and the Gaussian widths 
are allowed to float as free parameters in the likelihood fit 
and are extracted together with the yields. This strategy 
reduces the systematic error due to possible inaccuracies 
of the AE and rriEs description in Monte Carlo simula- 
tions. 

The ff)kd component is represented by a phenomeno- 
logical function with eight fixed parameters, all estimated 
from the distribution of S for the events in the mss side- 
band (Fig. I). 
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FIG. 5: The S distribution and fit for the events in the tties 
sideband of the data sample. 



The hbkd component is represented by the sum of an 
ARGUS and a Gaussian function, with parameters esti- 
mated from the distribution of tties for the events in the 
AEk and Ai?^ sidebands. 

The g components are each represented by a phe- 
nomenological function with seven fixed parameters. The 
parameters are estimated with Monte Carlo simulations 
for the TT and K hypotheses, and with events in the ra-^s 
sideband for the background case. A comparison of the D 
distributions in the three hypotheses shows that this vari- 
able, introduced by our procedure for factorizing PDFs, 
provides little discriminating power. 

From the maximum likelihood fit to the selected sam- 
ple we obtain iV^ = 52 ± 10, Nk = 1284 ± 37, and 
Nbkd = 819 ± 31. The correlation coefficient between TV^r 
and Nk is —0.04. The confidence level of the fit, defined 



as the probability to obtain a maximum value of the like- 
lihood smaller than the observed value, is 54%, estimated 
by Monte Carlo techniques. The statistical significance 
of the J/ip TT* signal, evaluated from the change in 

the maximum value of InL when we constrain = 0, 
is 7.0(7. 
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FIG. 6: The ln{P.^/PK) distribution for events in the data 
sample (after the subtraction of the background component 
in each bin) and from Monte Carlo simulations of 
J/'ipTv^{K^) events; the distributions are normalized to the 
yields extracted from the maximum likelihood fit. 



The distribution of ln(P^/PR-) for the sample, after 
subtraction of the background component in each bin, 
is shown in Fig. ^. The background distribution is nor- 
malized to the number of background events from the 
fit. The distribution of In^Pj^ / Pk) for simulated signal 
samples, normalized to the yields extracted from the like- 
lihood fit, is also shown. The distribution in Ai?^ for the 
events in the data sample with toes > 5.27GeV/c^ is 
shown in Fig. |^, along with the likelihood fit result. 

Possible biases in the fitting procedure were inves- 
tigated by performing the fit on simulated samples of 
known composition and of the same size as the data. The 
differences, A^ and Ak, between the extracted and the 
input values are consistent with 0. However we correct 
the yields for the observed deviations A^ = 1.1 ± 2.2 and 
Ak — —11.3 ± 8.8. The corrected yields are 51 ± 10 and 
1296 ±38 for J/ipTr"^ and J/i/jK^, respectively. 

The use of particle identification for the charged 
hadron has been investigated by adding to the likeli- 
hood, as an additional argument, the Cherenkov angle 
measured in the DIRC for this track. The PDFs for the 
variable 9c are determined from data and parameterized 
as Gaussian functions, with mean values and widths that 
depend on the momentum of the track. A fit with a mod- 
ified likelihood function is performed with the subsample 
of events where the particle identification information is 
available. The ratio of branching fractions is determined 
separately for the J/-0 (^"'"^~)/i^ and J/ip{e~^e~)h^ sam- 
ples. A detailed comparison, reported in Table |, shows 
that the addition of particle identification does not sig- 
nificantly change the statistical precision of the results, 
which are consistent to within l.Gcr. 
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FIG. 7: The AEt^ distribution for events with rriES > 
5.27GeV/c^ compared with the fit result (solid curve). 
The dotted curve represents the fitted contribution from 
the background alone, while the dashed curve repre- 
sents the fitted contributions from the sum of back- 
ground and J/ifj components. The PDFs of the AJ5-^ 
variable in the J/ip and background hypotheses have 
been obtained with a numerical integration of the Pi 
PDFs: pk{AE^) = J fK{x)gK{x- AE^)dx, pbkd{AE^) = 
/ fbkd{x + AEn)gtkd{x - AE„) dx. 



Based on the fitted event yields, we find the ratio of 
branching fractions to be 



= [3.91 ± 0.78{stat.) ± 0.19{syst. 



The dominant systematic error (0.17%) comes from the 
uncertainty in the correction factors, Ajr and Ak, due 
to the limited statistics of the simulated samples. The 
uncertainty in the fixed parameters of the PDFs, deter- 
mined by fits to simulated or non-signal data sets, affects 
several aspects of the likelihood fit: the characterization 
of the S and D distributions; the characterization of the 
Toes distribution for the background (including the frac- 
tion of peaking background events); and the fraction of 
signal events in the tails of the AE distribution. This 
uncertainty contributes 0.07% to the systematic error. 
Contributions due to any possible difference in the re- 
construction efficiencies for J/ipir^ and J/ipK^ events 
are found to be negligible, as are uncertainties due to 
inaccuracies in the description of the tails of the AE res- 
olution function. 

Our determination of the ratio of branching fractions 
is consistent with the expectation reported in |^ and 
with previous measurements ^, but has a substan- 
tially lower uncertainty than the world average value of 
(5.1 ±1.4)% §. 

To study direct CP-violation in these channels, we 



TABLE I: Measurements of B(B^ -> J/t/- 7r±)/B(B± 
J/tP K^) obtained with the original (fit 1) and a modified 
likelihood function (fit 2) that includes particle identification 
for h^. The error on the difference A between the two mea- 
surements is estimated as cta = \f\(j\ — a\\. 



sample 




fit 1 fit 2 


A/aA 


jj-i) (/i+/i 




(4.2 ± 1.0)% (4.7 ± 1.1)% 


1.1 


J/V'(e+e- 




(3.5 ±1.2)% (4.1 ±1.3)% 


1.2 



modify the likelihood function in Eq. g as follows: 

M 

L' = e-^'^^l[Y.P: {AEi , , TOEs^' , q')N, , (6) 

j=l i 

where q is the charge of . We factorize the PDFs as 
Pl{'^E„,P,TTT'ES,q) ^ Pi{AE^,p,mBs)ci{q) , (7) 

where Ci{q) is the probability for the final state charged 
hadron, in a certain hypothesis, to have charge q. The 
Ci can be written in terms of the CP-violating charge 
asymmetries Ai as 



c.(9) = ^[(l 



A)/+(g)±(i±A)/-(g)], (8) 



where 



riq) = 



N+ 



1 if g = ±1 
if g = -1, 

if q = ±1 

1 if q = -1. 



(9) 
(10) 
(11) 



The asymmetry observables Ai are allowed to float as 
free parameters in the likelihood fit and are extracted 
together with the yields. 

We impose additional requirements on the charged 
track in the events to be used in the fit, selecting only 
those tracks for which the tracking efficiency has been 
accurately measured from data. Tracks are required to 
have a polar angle in the range [0.41, 2.54] radians, to in- 
clude at least 12 DCH hits, to have pt > lOOMeV/c, and 
to point back to the nominal interaction point within 
1.5 cm in the vertical plane and within 3 cm along the 
longitudinal direction. The selected sample contains 982 
J/iph- and 970 B+ J/'tph+ candidates. 

From the maximum likelihood fit to the data sample we 
obtain A^ = 0.01 ± 0.22, Ak = -0.001 ± 0.030, Atkd = 
0.018±0.039. The correlation coefficient between Att and 
Ak is -0.03. 

The uncertainty in the fixed parameters of the PDFs, 
determined by fits to simulated or non-signal data sets, 
contributes 0.0056 and 0.0002 to the systematic error on 
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A-n and AK^ respectively. The difference in tracking effi- 
ciency between positively and negatively charged tracks 
- primarily pions - has been studied in hadronic events 
by comparing the independent SVT and DCH tracking 
systems. The corrections to the asymmetries A-n and 
Ak are negligible. The uncertainty on the corrections 
contributes 0.0026 and 0.0020 to the systematic error on 
A-^ and Ak, respectively. The fake asymmetry due to 
the different probability of interaction of and in 
the detector material before the DCH is estimated to be 
—0.0039. We correct Ak for this quantity and conserva- 
tively assume a contribution of 0.0039 to the systematic 
uncertainty. This represents the dominant systematic er- 
ror on Ak- a more careful evaluation of the materials 
and of jK^ cross-section differences will make it pos- 
sible to substantially reduce this contribution. 

We determine the CP-violating charge asymmetries to 

be 

A-, ^ 0.01±0.22(siai.)±0.01(syst.), 
Ak = 0.003 ±0.030(siai.)±0.004(s2;si.). 

These results are consistent with Standard Model expec- 
tations and with the measurement reported in . 

As a crosscheck, Ak has been determined also with a 
simple analysis based on the counting of J/ip 
signal events in the toes peak. The result is compat- 
ible with the likelihood fit analysis: Ak ~ 0.005 ± 
0.030(stat.)±0.004(syst.). 

We observe no evidence for CP-violation in 
J/ipn^ or B^ J/ifjK^ decays. These results are sta- 
tistically limited and can be expected to improve with 



additional data. 
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